AFRL-VA-WP-TP-2007-303 

FAULT  TOLERANT  OPTIMAL 
TRAJECTORY  GENERATION  FOR 
REUSABLE  LAUNCH  VEHICLES 
(PREPRINT) 


Patrick  J.  Shaffer,  I.  Michael  Ross,  Michael  W.  Oppenheimer, 
David  B.  Doman,  and  Kevin  P.  Bollino 


DECEMBER  2006 


Approved  for  public  release;  distribution  unlimited. 


STINFO  COPY 

This  is  a  work  of  the  U.S.  Government  and  is  not  subject  to  copyright  protection  in  the  United 
States. 


AIR  VEHICLES  DIRECTORATE 
AIR  FORCE  MATERIEL  COMMAND 
AIR  FORCE  RESEARCH  LABORATORY 
WRIGHT-PATTERSON  AIR  FORCE  BASE,  OH  45433-7542 


NOTICE  AND  SIGNATURE  PAGE 


Using  Government  drawings,  specifications,  or  other  data  included  in  this  document  for 
any  purpose  other  than  Government  procurement  does  not  in  any  way  obligate  the  U.S. 
Government.  The  fact  that  the  Government  formulated  or  supplied  the  drawings, 
specifications,  or  other  data  does  not  license  the  holder  or  any  other  person  or  corporation; 
or  convey  any  rights  or  permission  to  manufacture,  use,  or  sell  any  patented  invention  that 
may  relate  to  them. 

This  report  was  cleared  for  public  release  by  the  Air  Force  Research  Laboratory  Wright  Site 
(AFRL/WS)  Public  Affairs  Office  and  is  available  to  the  general  public,  including  foreign 
nationals.  Copies  may  be  obtained  from  the  Defense  Technical  Information  Center  (DTIC) 
(http://www.dtic.mil). 


AFRL-VA-WP-TP-2007-303  HAS  BEEN  REVIEWED  AND  IS  APPROVED  FOR  PUBLICATION  IN 
ACCORDANCE  WITH  ASSIGNED  DISTRIBUTION  STATEMENT. 


*//Signature// 

Michael  W.  Oppenheimer 
Electronics  Engineer 
Control  Design  and  Analysis  Branch 
Air  Force  Research  Laboratory 
Air  Vehicles  Directorate 


//Signature// 

Deborah  S.  Grismer 
Chief 

Control  Design  and  Analysis  Branch 
Air  Force  Research  Laboratory 
Air  Vehicles  Directorate 


//Signature// 

JEFFREY  C.  TROMP 
Senior  Technical  Advisor 
Control  Sciences  Division 
Air  Vehicles  Directorate 


This  report  is  published  in  the  interest  of  scientific  and  technical  information  exchange,  and  its 
publication  does  not  constitute  the  Government’s  approval  or  disapproval  of  its  ideas  or  findings. 

*Disseminated  copies  will  show  “//Signature//”  stamped  or  typed  above  the  signature  blocks. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


1 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  searching  existing  data 
sources,  gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of 
information,  including  suggestions  for  reducing  this  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports  (0704-0188),  1215  Jefferson  Davis 
Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  any  penalty  for  failing  to  comply  with  a  collection  of 
information  if  it  does  not  display  a  currently  valid  OMB  control  number.  PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 

1.  REPORT  DATE  (DD-MM-YY) 

December  2006 

2.  REPORT  TYPE  3.  DATES  COVERED  (From  -  To) 

Journal  Article  Preprint  06/0 1  /2004  -  12/1 8/2006 

4.  TITLE  AND  SUBTITLE 

FAULT  TOLERANT  OPTIMAL  TRAJECTORY  GENERATION  FOR  REUSABLE 
LAUNCH  VEHICLES  (PREPRINT) 

5a.  CONTRACT  NUMBER 

In-house 

5b.  GRANT  NUMBER 

5c.  PROGRAM  ELEMENT  NUMBER 

6220 IF 

6.  AUTHOR(S) 

Patrick  J.  Shaffer,  I.  Michael  Ross,  and  Kevin  P.  Bollino  (Naval  Postgraduate  School) 
Michael  W.  Oppenheimer  and  David  B.  Doman  (AFRL/VACA) 

5d.  PROJECT  NUMBER 

A03G 

5e.  TASK  NUMBER 

5f.  WORK  UNIT  NUMBER 

0B 

7.  PERFORMING  ORGANIZATION1 

Naval  Postgraduate  School 
Monterey,  CA  93943 

J  NAME(S)  AND  ADDRESS(ES) 

Control  Design  and  Analysis  Branch  (AFRL/VACA) 

Control  Sciences  Division 

Air  Vehicles  Directorate 

Air  Force  Materiel  Command,  Air  Force  Research  Laboratory 
Wright-Patterson  Air  Force  Base,  OH  45433-7542 

8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-VA-WP-TP-2007-303 

9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Vehicles  Directorate 

Air  Force  Research  Laboratory 

Air  Force  Materiel  Command 

Wright-Patterson  Air  Force  Base,  OH  45433-7542 

10.  SPONSORING/MONITORING 
AGENCY  ACRONYM(S) 

AFRL-VA-WP 

11.  SPONSORING/MONITORING 
AGENCY  REPORT  NUMBER(S) 

AFRL-VA-WP-TP-2007-303 

12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

Approved  for  public  release;  distribution  unlimited. 

13.  SUPPLEMENTARY  NOTES 

Journal  article  submitted  to  Guidance,  Control  and  Dynamics,  published  by  AIAA. 

This  is  a  work  of  the  U.S.  Government  and  is  not  subject  to  copyright  protection  in  the  United  States.  PAO  Case  Number: 
AFRL/WS  07-0100  (cleared  January  17,  2007).  Paper  contains  color. 

14.  ABSTRACT 

Reconfigurable  inner-loop  control  laws  improve  the  fault  tolerance  of  a  vehicle  to  control  effector  failures;  however,  in 
order  to  preserve  stability,  the  unfailed  effectors  may  be  deployed  to  off-nominal  positions  to  compensate  for  undesirable 
perturbations  caused  by  the  failed  effectors.  The  effectors  acting  under  the  influence  of  a  reconfigurable  control  law  can 
produce  significant  perturbations  to  the  nominal  forces  produced  by  the  wing  and  body  and  can  also  affect  the  range  of 
flight  conditions  over  which  the  vehicle  can  be  controlled.  Three  degree-of-freedom  (3  DOF)  dynamical  models  used  in 
trajectory  optimization  for  aerospace  vehicles  typically  include  wing-body  aerodynamic  force  effects  but  ignore  the 
aerodynamic  forces  produced  by  the  control  surfaces.  In  this  work,  a  method  for  including  these  trim  effects  as  well  as 
control  induced  trajectory  constraints  in  a  3  DOF  model  is  presented. 

15.  SUBJECT  TERMS 

Reconfigurable  Control,  Nonlinear  Control  Allocation,  Trajectory  Reshaping 

16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION 

18.  NUMBER  OF 

19a.  NAME  OF  RESPONSIBLE  PERSON  (Monitor) 

a.  REPORT 

Unclassified 

b.  ABSTRACT 

Unclassified 

c.  THIS  PAGE 

Unclassified 

OF  ABSTRACT: 

SAR 

PAGES 

42 

Michael  W.  Oppenheimer 

19b.  TELEPHONE  NUMBER  (Include  Area  Code) 

N/A 

Standard  Form  298  (Rev.  8-98) 

Prescribed  by  ANSI  Std.  Z39-18 


1 


Fault  Tolerant  Optimal  Trajectory  Generation  for  Reusable 


Launch  Vehicles 

Patrick  J.  Shaffer  * 

I.  Michael  Ross  t 
Michael  W.  Oppenheimer  * 

David  B.  Doman  § 

Kevin  P.  Bollino  ' 

Reconfigurable  inner-loop  control  laws  improve  the  fault  tolerance  of  a  vehicle  to  control 
effector  failures;  however,  in  order  to  preserve  stability,  the  unfailed  effectors  may  be  de¬ 
ployed  to  off-nominal  positions  to  compensate  for  undesirable  perturbations  caused  by  the 
failed  effectors.  The  effectors  acting  under  the  influence  of  a  reconfigurable  control  law  can 
produce  significant  perturbations  to  the  nominal  forces  produced  by  the  wing  and  body  and 
can  also  affect  the  range  of  flight  conditions  over  which  the  vehicle  can  be  controlled.  Three- 
degree-of- freedom  (3  DOF)  dynamical  models  used  in  trajectory  optimization  for  aerospace 
vehicles  typically  include  wing-body  aerodynamic  force  effects  but  ignore  the  aerodynamic 
forces  produced  by  the  control  surfaces.  These  so-called  trim  effects  are  a  normal  com¬ 
ponent  of  six-degree-of-freedom  models;  however,  such  models  have  been  considered  to  be 
excessively  cumbersome  for  the  purposes  of  trajectory  optimization.  In  this  work,  a  method 
for  including  these  trim  effects  as  well  as  control  induced  trajectory  constraints  in  a  3  DOF 
model  is  presented.  The  method  makes  use  of  nonlinear  control  allocation  methods  to  deter¬ 
mine  the  control  effector  positions  required  to  rotationally  balance  the  vehicle  at  each  point 
in  the  flight  envelope  and  uses  that  information  to  compute  force  perturbations  and  new 
control  induced  trajectory  constraints.  The  method  estimates  the  steady  state  response  of  a 
reconfigurable  inner-loop  control  law  by  solving  for  the  control  surface  deflections  required 
to  achieve  rotational  equilibrium  over  a  range  of  flight  conditions.  This  is  accomplished  by 
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passing  knowledge  of  surface  failures  to  a  nonlinear  control  allocation  algorithm  in  order  to 
estimate  the  control  surface  deflections  required  to  balance  the  vehicle.  Upon  solving  for  the 
required  control  surface  deflections,  the  contributions  of  the  control  effectors  to  the  vehicle 
lift  and  drag  are  computed,  as  well  as  3  DOF  compatible  constraints  that  define  regions  of 
the  flight  envelope  where  the  vehicle  does  not  have  sufficient  control  power  to  rotationally 
balance  the  vehicle.  An  example  is  presented  that  includes  control  failure  effects  in  a  max¬ 
imum  downrange  trajectory  optimization  problem  for  the  flight  of  an  unpowered  reusable 
launch  vehicle  passing  through  the  atmosphere.  Constrained  trajectories  are  planned  by  solv¬ 
ing  the  resulting  optimal  control  problem  using  a  Legendre  pseudospectral  method.  State 
dependent  constraints  resulting  from  a  control  failure  are  included  in  the  nonlinear  optimal 
control  problem.  With  accurate  knowledge  of  the  available  maneuverability  envelope  and 
aerodynamics  of  the  failed  vehicle,  it  is  demonstrated  that  the  methods  hold  the  potential 
to  significantly  enhance  safety  margins  even  when  control  effectors  have  failed. 

Introduction 

Autonomous  reusable  launch  vehicles  (RLV)  are  being  pursued  as  low-cost  alternatives 
to  expendable  launch  vehicles  and  the  Shuttle.  The  employment  of  autonomous,  reusable 
launch  vehicles  requires  additional  robustness  in  guidance  and  control  to  fulfill  the  role  of 
an  adaptive  human  pilot  in  the  event  of  failures  or  unanticipated  conditions.  Hence,  new 
guidance  strategies  that  are  able  to  identify  and  adapt  to  vehicle  failures  during  the  flight 
and  still  safely  return  the  vehicle  to  earth  are  required.1  New  advanced  guidance  and  control 
algorithms  (AG&C)  have  been  developed  that  adapt  and  attempt  to  recover  vehicles  that 
have  suffered  failures  or  damage.  A  survey  of  historical  launch  vehicle  failures  that  resulted 
in  vehicle  loss  or  degraded  mission  capabilities  concluded  that  AG&C  methods  could  improve 
mission  performance  or  recover  a  reusable  launch  vehicle  that  experienced  a  similar  failure.2 

In  this  work,  the  framework  for  computing  optimal  trajectories  for  a  nominal  and  failed 
vehicle  is  developed.  The  method  presented  here  is  designed  to  be  used  in  conjunction  with 
reconfigurable  inner-loop  flight  control  laws.  While  such  control  methods  improve  the  fault 
tolerance  of  a  vehicle  to  control  effector  failures;  the  unfailed  effectors  may  be  deployed  in 
unusual  ways  to  compensate  for  undesirable  perturbations  caused  by  the  failed  effectors. 
Reconfigurable  control  laws  can  therefore  produce  significant  perturbations  to  the  nominal 
forces  produced  by  the  wing  and  body  and  can  also  affect  the  range  of  flight  conditions  over 
which  the  vehicle  can  be  controlled.  Three-degree-of-freedom  (3  DOF)  dynamical  models 
used  in  trajectory  optimization  for  aerospace  vehicles  typically  include  wing-body  aerody¬ 
namic  force  effects  but  ignore  the  aerodynamic  forces  produced  by  the  control  surfaces.  These 
so-called  trim  effects  are  a  normal  component  of  six-degree-of-freedom  models;  however,  such 
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models  have  historically  been  considered  less  attractive  than  3  DOF  models  for  the  purposes 
of  trajectory  optimization.  A  method  for  including  trim  effects  as  well  as  control  induced 
trajectory  constraints  in  a  3  DOF  model  is  presented.  The  method  makes  use  of  recent 
advances  in  nonlinear  control  allocation  methods  to  determine  the  control  effector  positions 
required  to  rotationally  balance  the  vehicle  at  each  point  in  the  flight  envelope  and  uses  that 
information  to  compute  force  perturbations  and  new  control  induced  trajectory  constraints. 
The  method  estimates  the  steady  state  response  of  a  reconfigurable  inner-loop  control  law 
by  solving  for  the  control  surface  deflections  required  to  achieve  rotational  equilibrium  over 
a  range  of  flight  conditions.  This  is  accomplished  by  passing  knowledge  of  surface  failures  to 
a  nonlinear  control  allocation  algorithm  in  order  to  estimate  the  control  surface  deflections 
required  to  balance  the  vehicle.  Upon  solving  for  the  required  control  surface  deflections,  the 
contributions  of  the  control  effectors  to  the  vehicle  lift  and  drag  are  computed,  as  well  as  3 
DOF  compatible  constraints  that  define  regions  of  the  flight  envelope  where  the  vehicle  does 
not  have  sufficient  control  power  to  rotationally  balance  the  vehicle.  Such  constraints  are 
state  dependent  and  are  utilized  in  an  optimal  control  problem  to  replan  trajectories  that 
avoid  regions  of  the  flight  envelope  where  the  vehicle  cannot  maintain  rotational  equilibrium. 

Trajectory  and  abort  planning  functions  use  a  variety  of  optimality  criteria  as  implicit 
requirements.  For  example,  in  abort  planning  for  unpowered  re-entry  vehicles,  the  size  of  the 
footprint  or  reachable  region  on  the  surface  of  the  earth  is  an  important  factor  to  consider. 
Maximizing  vehicle  footprint  under  failure  conditions  is  important  in  order  to  allow  one  to 
select  from  as  many  landing  sites  as  possible  during  exigency  operations.  Other  trajectory 
generation  and  guidance  problems  include  ascent,  re-entry,  terminal  area  energy  management 
and  approach  and  landing.  Mease,  et.al.,3  developed  a  three  dimensional  entry  trajectory 
planning  method.  Shen  and  Lu4  develop  a  methodology  to  design  3  DOF  entry  trajectories 
and  also  describe  a  new  automated  lateral  guidance  logic  based  on  crossrange.5  Some  of 
these  methods  focus  on  designing  reference  trajectories  on-board  in  order  to  improve  perfor¬ 
mance. 3,4>6> 7  While  these  studies  did  not  directly  account  for  the  effects  of  control  effector 
failures  and  reconfigurable  control  system  response,  these  methods  could  conceivably  make 
use  of  the  methods  developed  in  this  work  that  update  the  3  DOF  aerodynamic  model  to 
include  trim  force  effects  and  flight  envelope  constraints.  Related  work  has  been  accom¬ 
plished  that  integrates  adaptive  guidance  and  control  methods  with  trajectory  modification 
methods  to  automatically  respond  to  failures8,23  and  design  abort  trajectories  online.9 

For  the  purpose  of  illustrating  the  method  of  incorporating  control  surface  failure  effects 
into  a  3  DOF  aerodynamic  model  and  estimating  control  failure  induced  flight  envelope 
constraints,  we  consider  the  problem  of  determining  the  maximum  downrange  capability  for 
an  unpowered  RLV  under  nominal  and  failed  conditions.  Under  failure  conditions,  the  safety 
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margins  of  traditional  feedback  control  laws  can  diminish  because  they  are  not  tuned  for  the 
failed  vehicle  or  because  they  are  based  on  inadequate  dynamical  models.10  Oversimplified 
dynamical  models  that  do  not  exploit  the  full  maneuverability  envelope  of  the  vehicle  can 
compromise  safety.11  By  incorporating  the  forces  produced  by  the  control  effectors  into 
the  3  DOF  models  used  for  trajectory  optimization,  one  can  reduce  aerodynamic  model 
inaccuracies  and  estimate  regions  in  the  flight  envelope  that  should  be  avoided  due  to  a  lack 
of  control  power  and  the  inability  to  maintain  rotational  equilibrium  at  critical  points  in 
the  flight  envelope.  Therefore,  determining  optimal  trajectories  under  failure  conditions  is 
critical  to  safety.  For  example,  a  conventional  footprint  generation  algorithm  may  incorrectly 
assess  the  feasibility  of  reaching  a  critical  landing  site,  while  a  more  sophisticated  algorithm 
could  generate  revised  optimal  trajectories  based  on  accurate  failure  effect  information  in 
order  to  safely  recover  a  vehicle  at  a  feasible  landing  site.12  Thus,  failure  to  perform  trajectory 
optimization  using  accurate  aerodynamic  models  and  flight  envelope  constraint  estimates  for 
the  failed  vehicle  could  lead  to  incorrect  conclusions  regarding  landing  site  feasibility  during 
online  abort  planning. 

The  class  of  failures  considered  in  this  work  are  control  surface  failures  whose  effects 
can  be  described  using  the  nominal  aerodynamic  database.  Nominal  trajectories  are  not 
necessarily  feasible  following  a  control  surface  failure.  The  challenge  is  to  reconfigure  the 
inner-loop  control  system  to  use  the  healthy  control  effectors  to  compensate  for  the  undesired 
effects  of  the  failed  control  effector  and  to  replan  a  new  trajectory  based  upon  aerodynamic 
data,  that  includes  off-nominal  trim  effects,  to  achieve  an  acceptable  end  condition.  Despite 
the  fact  that  this  is  an  extremely  difficult  problem,  recent  research  on  pseudospectral  (PS) 
methods  for  optimal  control  indicates  that  solving  such  problems  are  indeed  within  reach.13, 14 
In  this  paper,  we  describe  preliminary  steps  in  achieving  this  ultimate  goal. 

This  work  builds  upon  the  recent  research  of  Fahroo  and  Doman12  in  demonstrating 
the  application  of  a  PS  method  to  determine  optimal  trajectories  for  an  RLV  for  nominal 
conditions  as  well  as  for  off-nominal  conditions  caused  by  vehicle  control  surface  failures.  To 
this  end,  this  paper  integrates  several  different  aspects  of  the  guidance  and  control  problem. 
X-33  vehicle  aerodynamic  data15  is  used  in  a  3  DOF  dynamical  model  to  generate  the 
trajectories.  The  main  thrust  of  this  work  is  the  inclusion  of  state-dependent  Mach  and 
angle  of  attack  constraints  in  the  optimization  problem.  These  constraints  are  computed 
by  forming  a  trim  deficiency  map,  which  displays  regions  in  the  Mach-angle  of  attack  space 
where  the  vehicle  can  be  rotationally  trimmed.  The  boundary  of  the  trimmable  region  is 
used  as  a  constraint  in  the  optimization  problem,  requiring  that  any  reshaped  trajectories 
he  in  the  trimmable  region.  The  example  utilized  here  is  the  problem  of  determining  the 
maximum  downrange  for  a  footprint  analysis.  Heating  and  dynamic  pressure  constraints 
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are  not  included  in  the  optimization  problem  (normal  force  is  included),  however,  these  can 
be  easily  incorporated  with  the  framework  presented  here.  In  essence,  the  work  presented 
here  is  general  and  can  accommodate  a  wide  range  of  trajectory  problems.  With  suitable 
modifications  to  the  optimization  problem  (changes  to  the  cost  and  constraints),  entry  and 
approach  and  landing  guidance  problems  and  heating  and  dynamic  pressure  constraints  can 
be  considered. 


Constructing  the  Dynamical  Model:  Motivation  and  Issues 

For  guidance  purposes,  it  is  common  practice  to  consider  a  3  DOF  model  that  describes 
the  point-mass  motion.  The  angle  of  attack,  a,  and  bank  angle,  a,  are  selected  to  be  the 
control  inputs.  A  3  DOF  model  for  flight  over  a  flat  earth  is  given  by:16 


x 

y 

z 

V 

7 

X 


V  cos  %  cos  7 

V  sin  x  cos  7 

V  sin  7 

Cd(M,  a)ArefpV2 

- - - g  sin  7 

Am 

Cl(M,  a) ArefpV  cos  a  g  cos  7 

2m  V 

Cl(M ,  a)ArefpV  sincr 

2  m 


(1) 

(2) 

(3) 

(4) 

(5) 

(6) 


where  x  (down  range),  y  (cross-range),  and  2  (altitude)  are  the  vehicle’s  position,  7  is  the 
flight  path  angle,  x  is  the  azimuth  angle,  a  is  the  bank  angle,  V  is  the  velocity  magnitude,  p 
is  the  freestream  mass  density,  Aref  is  a  characteristic  area  for  the  body,  and  Cl  and  Cd  are 
the  lift  and  drag  coefficient  functions  that  depend  upon  the  Mach  number,  M ,  and  angle  of 
attack,  a. 

In  order  to  ensure  smoothness  of  the  controls  q-(-)  and  <r(-),  we  define  the  rates  d(-)  and 
&(•)  as  the  control  inputs  and  add  two  additional  states  to  the  3  DOF  equations  of  motion: 


a{t)  =  ua(t )  (7) 

a(t)  =  ua(t )  (8) 


where  u  :=  (ua,ua)  G  M2  are  the  controls  that  must  be  selected  from  some  compact  set, 
UcR2.  In  this  work,  we  use 


U  :=  {(ua,uc r)  :  \ua\  <  40  deg/s ,  \ua\  <  40  deg/s} 


(9) 
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It  is  worth  noting  that  in  the  absence  of  this  specification,  a(-)  and  a(-)  may  exhibit  highly 
dynamic  behavior  depending  upon  the  type  of  guidance  problem.15,17  Thus,  by  adding  these 
rate  dynamics  to  the  3  DOF  model,  we  arrive  at  an  eight  dimensional  state  vector,  that  is, 
x  e  M8. 

The  vehicle  lift  and  drag  coefficients,  Cl  and  C-o ,  are  functions  of  Mach  number  and 
angle  of  attack  and  can  be  broken  down  into  a  component  from  the  base  or  wing-body  and 
increments  due  to  the  deflection  of  aerodynamic  control  surfaces:15 


CL  (M,  a)  =  CLo  (M,  a)  +  CLg, .  (M,  a,  S*(M,  a)) 
CD  (M,  a)  =  CDo  (M,  a)  +  CDg,  (M, «,  6*(M,  a)) 


where  ( M,a )  i— >  ( Cl0,Cd0 )  are  based  upon  static  table-look  up  data  while  (M,a)  i— > 
( Cls »,  Cds*)  are  based  upon  an  algorithmic  map,  (M,  a)  i— >  5*,  and  <5*  is  a  solution  to  a  con¬ 
trol  allocation  algorithm  used  to  trim  the  vehicle.  The  quantities,  Cl0  (M,  a)  and  Coa  (M,  a) 
are  base  lift  and  drag  coefficients,  while  Cls*  (M,  a,  8*(M,  a))  and  Cds*  (M,  a,  8*(M,  a))  are 
the  control  induced  lift  and  drag  increments.  The  quantity  6  =  6*  corresponds  to  a  trim 
vector  of  inner-loop  control  variables  that  establish  rotational  equilibrium  at  a  given  flight 
condition.  Nominally,  there  exist  8  aerodynamic  control  surfaces  on  the  X-33:  left  and  right 
rudders,  body  flaps,  inboard  and  outboard  elevons.  The  X-33  is  shown  in  Fig.  1.  These 
control  surface  deflection  variables,  collectively  known  as  the  effector  vector,  are  given  by 


<5 


ElevoriRj,  ElevoriR0,  ElevoriLn  ElevoriL0,  FlapR,  FlapL ,  RudderR, 


Rudder  l 


(11) 


In  Eq.  11,  the  subscripts  are  defined  as  i?/=right  inboard,  Rq =right  outboard,  L/=left 
inboard,  Lo=left  outboard,  i?=right,  and  L= left.  The  trim  control  deflection  vector  d*  is 
obtained  by  solving  a  nonlinear  control  allocation  problem  that  enforces  actuator  limits  and 
accommodates  locked  or  floating  aerodynamic  surface  failures. 

It  should  be  noted  that  in  this  exposition,  the  effect  of  sideslip  angle  on  the  lift  and  drag 
forces  is  ignored  for  the  sake  of  clarity.  Adding  a  third  independent  variable  would  result  in 
four-dimensional  trim  constraint  regions,  which  would  be  difficult  to  view  and  analyze.  This 
manuscript  considers  failures  for  which  there  exists  sufficient  control  power  to  maintain  zero 
sideslip  angle.  Nevertheless,  the  methods  presented  here  are  general  in  that  flight  at  non-zero 
sideslip  conditions  can  be  accommodated;  however,  the  lift,  drag,  and  moment  coefficients 
become  functions  of  Mach,  angle  of  attack,  and  sideslip  thereby  making  visualization  of  trim 
maps  difficult. 
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Control  Allocation  in  Failure  Modes 

The  RLV  control  surfaces  are  position  limited  and  have  maximum  and  minimum  deflec¬ 
tions.  In  this  work,  6  is  limited  by 


d  m.in. 


-30, 


dr,,  n  ,• 


-30,  -30,  -30,  -15,  -15,  -60, 
30,  30,  30,  30,  26,  26,  30,  60 


(12) 


where  the  numbers  are  in  degrees.  Numerous  reconfigurable  inner-loop  control  methods  exist 
whose  purpose  is  to  maintain  control  of  a  vehicle  (when  physically  possible)  in  the  event  of 
a  failure  or  damage  by  using  the  healthy  control  surfaces  to  compensate  for  the  effects  of 
the  failed  surfaces. 1’8, 18-22  Under  some  failure  conditions,  reconfigurable  inner-loop  control 
is  insufficient  to  recover  a  vehicle  that  also  relies  on  an  autonomous  guidance  system  that 
follows  a  prescribed  trajectory.23  Significant  changes  to  the  aerodynamics  and  a  possible  lack 
of  control  power  may  force  one  to  determine  a  new  trajectory  based  on  knowledge  of  the  effect 
of  a  failure.  Failures  may  result  in  changes  in  the  constraints  on  the  state,  cost,  control,  and 
path  arguments  that  correspond  to  the  problem  of  computing  a  new  optimal  trajectory  for  a 
failed  vehicle.  A  replanned  trajectory  must  be  computed  that  provides  the  ability  to  complete 
a  mission  to  a  terminal  objective  in  order  to  safely  recover  the  vehicle.  Traditional  methods 
store  the  trajectory  data  for  a  limited  set  of  pre-planned  engine-out  or  control  surface  failure 
scenarios  and  therefore  require  significant  pre-planning,  forethought,  and  data  storage  for 
implementation.  While  data  storage  problems  are  beginning  to  disappear,  conceiving  all 
possible  failure  combinations  and  designing  trajectories  for  each  case  is  at  best,  a  formidable 
and  painstaking  task.  An  appealing  alternative  is  to  develop  a  fault  tolerant  autonomous 
guidance,  control,  and  trajectory  shaping  system  that  can  be  executed  in-flight  and  respond 
to  failures  as  they  occur.  In  support  of  this  goal,  we  describe  the  development  of  an  online 
algorithm  that  can  provide  a  trajectory  reshaping  algorithm  with  an  updated  aerodynamic 
database  that  accurately  represents  the  effects  of  locked  or  floating  aerodynamic  surfaces, 
and  provides  a  constraint  map  that  describes  the  feasible  range  of  Mach  number  and  angle  of 
attack  for  which  the  vehicle  can  be  rotationally  trimmed.  Thus,  the  algorithm  allows  one  to 
include  6  DOF  trim  effects  and  constraints  in  the  3  DOF  dynamical  model  described  earlier. 

It  is  assumed  that  the  vehicle  has  sufficient  control  authority  to  trim  in  symmetric  flight 
(zero  sideslip  angle  )  and  that  the  total  lateral  wing-body  force  and  moment  coefficients  can 
be  driven  to  zero.  We  also  assume  that  the  body- axis  angular  velocity  vector  is  small  such 
that  damping  derivative  contributions  to  the  vehicle  moment  are  insignificant.  Thus,  the 
base  moments  for  this  study  are  functions  of  Mach  number  and  angle  of  attack  only.  The 
necessary  conditions  to  rotationally  trim  the  vehicle  are  that  the  moments  resulting  from  all 
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control  surface  deflections  (failed  and  healthy)  must  be  equal  and  opposite  to  the  moments 
produced  by  the  wing-body.  Thus, 


(  Crms(M,a,5)  \ 
Cms(M,a,6) 
Cyms(M,a,6 )  ) 


0  \ 

-Cmo(M,  a) 

0  ) 


(13) 


where  Crms(M,a,6),  Cm6(M,  a,  6),  and  Cyms(M,a,6)  are  the  rolling,  pitching,  and  yaw¬ 
ing  moment  coefficients  produced  by  the  control  effectors  and  Cnio  (M,  a)  is  the  wing-body 
pitching  moment. 

A  nonlinear  control  allocation  algorithm24  is  used  to  determine  the  effector  displacement 
vector  that  satisfies  Eq.  13  and  attempts  to  simultaneously  minimize  the  1-norm  of  the 
displacement  vector.  This  vector  is  nominally  subject  to  the  limits  of  Eq.  12  or  by  a  modified 
set  of  position  limits  that  reflect  the  current  health  of  the  effector  suite.  The  control  allocation 
problem  uses  an  iterative  algorithm25  that  queries  the  full  6  DOF  aerodynamic  database. 
The  constrained  control  allocation  problem  is  given  by: 


min  Jd 
5 


=  min 
5 


0 

C'rmsj^i  (Mj ,  OL\,  ^  j,i) 

(Mj  i  (yi ) 

— 

Crr iSj^i  {Mj ,  OLi,  $j,i) 

0 

_  C'ymsj^  {Mj ,  Sj,i)  _ 

subject  to 


<  8  <  8 


max 


(14) 


(15) 


at  each  Mach- a  point  in  the  aerodynamic  database.  Notice  that  the  upper  and  lower  bounds 
on  6  are  dependent  on  the  operating  mode  of  the  vehicle.  A  vehicle  operating  under  nominal 
conditions  has  the  full  range  of  deflections  available,  while  a  vehicle  operating  with  a  control 
surface  failure  has  a  subset  of  the  nominal  displacement  vector’s  range.  A  nonzero  value  for 
Eq.  14  means  that  there  is  insufficient  control  power  to  rotationally  trim  the  vehicle  (it  is 
trim  deficient),  and  that  the  vehicle  would  be  at  risk  of  departing  controlled  flight  at  that 
particular  flight  condition.  The  trim  deficient  regions  in  Mach- a-  space  can  be  mapped  and 
identified  as  flight  envelope  boundaries  that  constrain  the  trajectory.  If  the  control  alloca¬ 
tion  algorithm  is  able  to  rotationally  trim  the  vehicle,  any  excess  control  power  can  be  used 
to  maneuver  the  vehicle.  The  control  allocator  employs  a  second  control  allocation  opti¬ 
mization  algorithm  to  minimize  the  deviation  from  a  minimum  control  deflection  condition. 
This  condition  is  represented  by  an  eight-element  preference  vector  set  to  zero  and  roughly 
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corresponds  to  a  minimum  drag  condition.  The  second  optimization  is 


mjn  Js  =  min  ||W  (S  -  <SP)  ||x  (16) 

0  0 

subject  to 

&min  —  ^  —  &max  (17) 

and 


Crms(M,  a 

/ 

0 

\ 

Cms(M,a , 

S) 

= 

-Cmo(M,a) 

\  Cyms(M,a. 

.*) ) 

V 

0 

J 

In  Eq.  16,  W  is  a  weighting  matrix  used  to  weight  the  importance  of  driving  an  effector  to 
its  preferred  location.  In  the  case  where  the  vehicle  is  operating  with  excess  control  power, 
the  effector  displacement  vector  is  the  output  of  the  second  optimization  routine,  otherwise 
S  minimizes  the  trim  deficiency.  Here  we  set  the  preference  vector  6P  =  0  since  this  roughly 
corresponds  to  minimizing  control  effector  drag  which  is  compatible  with  the  higher  level 
objective  of  maximizing  downrange.  The  effector  displacement  vector  that  is  the  solution  to 
the  control  allocation  problem  is  represented  as  (5*  to  delineate  it  from  the  possibly  multiple 
combinations  of  control  surface  deflections  that  are  available  to  achieve  the  same  desired 
control  but  are  not  optimized  for  minimum  control  deflection. 

It  is  important  that  a  nonlinear  control  allocation  algorithm  be  used  to  solve  for  the 
deflection  vector  because  the  moments  are  nonlinear  functions  of  control  surface  deflection. 
While  linear  methods  may  be  suitable  for  inner-loop  flight  control  applications,  they  fail  to 
capture  the  global  character  of  nonlinear  moment-deflection  curves  and  may  lead  to  erroneous 
trim  deflection  vectors  that  fail  to  trim  the  vehicle.  The  piecewise  linear  control  allocation 
method  of  Bolender24  was  used  in  this  case  and  has  been  found  to  yield  acceptable  results 
for  cases  involving  separable  nonlinearities. 

Incorporating  Trim  Deficiency  Maps  In  Trajectory  Replanning 

The  vehicle  is  assumed  to  operate  in  the  nominal  condition  until  a  failure  has  been 
identified  through  an  online  fault  detection  or  identification  algorithm,  at  which  point  a  new 
lift  and  drag  table  corresponding  to  the  failed  vehicle  configuration  is  generated.  Once  a 
failure  has  been  identified,  the  control  allocator  generates  a  new  set  of  control  deflections 
for  the  failed  vehicle  and  uses  the  aerodynamic  database  to  determine  the  new  trim  lift 
and  drag  contributions.  In  the  process  of  computing  the  control  deflections  for  the  failed 
vehicle,  points  in  the  Mach-a  envelope  for  which  the  vehicle  cannot  be  trimmed  are  identified. 
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These  trim  deficiency  regions  must  be  avoided  and  are  thus  used  as  path  constraints  for  the 
trajectory  replanning.  A  point  in  the  Mach- a  space  is  called  trim  deficient  if  no  control 
effector  vector,  8,  can  be  found  to  identically  satisfy  Eq.  13.  The  aerodynamic  data  table  is 
then  changed  to  account  for  the  current  failed  condition  and  a  new  trajectory  is  generated 
for  the  new  aerodynamics  and  path  constraints.  Figures  2  and  3  show  the  trim  deficiency 
maps  for  a  nominal  case  and  a  case  of  both  bodyflaps  failed  at  26°.  In  three-dimensions, 
Fig.  3  becomes  the  plot  shown  in  Fig.  4.  Notice,  from  Fig.  2,  that  the  maximum  trim 
deficiency  value  is  on  the  order  of  machine  precision  so  that  the  entire  range  of  combinations 
of  angles  of  attack  and  Mach  numbers  is  considered  to  be  available  in  the  nominal  (unfailed) 
case.  Hence,  the  unfailed  vehicle  can  be  trimmed  over  the  entire  Mach- a  space.  When  the 
bodyflaps  are  locked  at  26°,  the  envelope  describing  allowable  Mach  numbers- angles  of  attack 
is  now  severely  reduced  (see  Figs.  3  and  4).  Points  in  the  Mach- a-  envelope  where  the  pitch 
deficiency  is  non-zero  indicate  that  there  is  insufficient  control  power  to  rotationally  balance 
the  vehicle  and  that  vehicle  control  would  be  in  jeopardy  at  such  conditions.  This  establishes 
a  restricted  region  that,  if  possible,  feasible  trajectories  should  avoid.  Ideally,  the  trajectory 
will  have  a  clear  path  in  Mach- a  space  from  the  initial  conditions  to  the  final  conditions.  A 
clear  path  is  defined  as  a  path  from  one  point  defined  by  the  flight  condition  (Mach-a)  to 
another  point  (Mach- a)  with  no  intervening  regions  of  unacceptable  trim  deficiency.  Should 
a  path  cross  into  or  through  a  region  of  pitch  deficiency,  the  trajectory  is  declared  infeasible 
because  the  vehicle  cannot  be  rotationally  balanced  over  part  of  the  candidate  trajectory. 
Failure  conditions  can  result  in  situations  where  no  clear  path  exists.  In  such  cases,  one 
might  select  a  path  where  the  least  severe  moment  deficiency  is  encountered. 

Having  obtained  the  control  effector  settings  that  rotationally  balance  the  vehicle,  the 
total  lift  and  drag  acting  on  the  vehicle  can  be  computed.  Utilizing  8*(M,  a)  from  the  control 
allocation  problem,  the  aerodynamic  database  can  be  queried  to  determine  the  control  surface 
induced  lift  and  drag  contributions.  Figures  5  and  6  show  the  nominal  and  failed  lift  and 
drag  maps  over  the  entire  Mach- a  range  of  interest.  In  this  case,  the  failure  of  the  bodyflaps 
at  26°  does  not  significantly  affect  the  lift  of  the  vehicle,  however,  large  changes  are  seen  in 
the  drag.  The  constraint  boundary  and  updated  lift  and  drag  maps  are  made  available  to  a 
trajectory  reshaping  algorithm  to  ensure  that  a  feasible  trajectory,  if  one  exists,  is  computed. 

The  algorithm  to  compute  pitch  deficiency,  lift,  and  drag  maps  and  to  determine  the 
range  of  trimmable  angle  of  attack  is  summarized  as  follows:25 

1.  Define  a  grid  for  Mach  and  Angle  of  Attack,  including  lower  and  upper  bounds  and 
step  size 

2.  Initialize  8 
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3.  Loop  1  — >  For  j  =  1  to  Number  of  Machs 

4.  Loop  2  — ►  For  i  =  1  to  Number  of  Angles  of  Attack 

5.  Compute  the  wing-body  pitching  moment  coefficient,  Cmoji(Mj,ati ) 

6.  Solve  a  control  allocation  problem  to  find  6*t  which  satisfies  Eq.  13 

ii  r  o 


7.  Compute  trim  deficiency  at  each  point,  i.e. ,  deficiency (j,i) 


(-'uioj.i  {Mj ,  cq) 

0 


Cq,  & j.i) 

(-  'msj.i  {Mj ,  Ot ,  8j  i) 

_  Cymsj,i(Mj,Oii,6jti)  _  2 

8.  If  deficiency (j,i)  =  0,  the  vehicle  can  be  trimmed  at  the  corresponding  Mach-O:  combi¬ 
nation 

9.  Compute  drag  and  lift  at  each  data  point  by  substituting  <5 A  into  the  aerodynamic 
database 

10.  Increment  Mach  and/or  Alpha 

11.  End  Angle  of  Attack  loop 

12.  End  Mach  loop 


This  algorithm  yields  the  control  deficiency  map  as  well  as  the  lift  and  drag  maps.  Each 
map  is  valid  for  all  operating  conditions  for  which  a  model  is  available. 

Note  that  the  previous  discussion  does  not  consider  any  other  limitations  such  as  normal 
force,  heating,  or  dynamic  pressure  on  the  trajectory.  These  constraints  may  further  reduce 
the  area  available  for  a  feasible  trajectory  and  can  be  taken  into  account  when  solving  the 
optimal  control  problem  by  appending  additional  inequality  constraints.  It  is  possible  that, 
in  some  cases,  these  additional  constraints  will  completely  eliminate  a  clear  path.  It  is 
easier  to  visualize  the  optimal  trajectory  with  constraints  as  a  region  through  Mach- a  space 
for  which  the  vehicle  can  be  rotationally  trimmed  without  violating  normal  force,  heating, 
or  dynamic  pressure  constraints.  If  a  failure  causes  this  corridor  to  close,  a  clear  path  to 
the  end  mission  may  not  exist.  Therefore,  feasible  trajectories  are  not  available  and  flight 
termination  may  be  the  only  option. 
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Trajectory  Replanning  Optimal  Control  Problems 

A  typical  trajectory  computation  problem  considers  a  vehicle,  at  some  initial  position, 
and  seeks  to  generate  an  optimal  trajectory  that  brings  it  to  a  final  state  specified  by  an 
endpoint  condition,  subject  to  constraints.  The  resulting  optimal  trajectory  minimizes  some 
performance  index,  which  can  include  some  function  of  the  vehicle  states  and/or  control 
inputs.  Typical  performance  indices  are  maximizing  downrange,  maximizing  cross  range, 
minimizing  control  effort,  or  some  weighted  combination  of  these  indices.  In  this  section,  we 
illustrate  how  to  incorporate  the  effects  of  control  failures,  using  the  methods  discussed  in 
the  previous  sections,  into  the  solution  of  a  maximum  downrange  problem  for  an  unpowered 
RLV. 

As  discussed  in  the  previous  sections,  a  3  DOF  dynamic  model,  that  incorporates  infor¬ 
mation  derived  from  a  6  DOF  model,  will  be  used.  In  particular,  we  incorporate  the  effects 
of  aerodynamic  surface  induced  lift  and  drag  from  a  high-fidelity  aerodynamic  database  and 
we  determine  the  region  in  Mach- a  space  over  which  rotational  equilibrium  can  be  achieved. 
In  order  to  determine  a  trajectory  that  maximizes  downrange,  we  select  the  following  cost 
function: 

J[x(-), u(-), to, */]  :=  -x{tf)  (19) 

where  the  elements  of  x  are  given  by, 


x  = 


xyzM^xa(J 


T 


(20) 


The  initial  conditions  for  the  example  mission  are 


T  -  \ 

(  0 

Vi 

0 

A 

125, 000ft 

Mi 

8 

7* 

-1.3° 

X* 

0 
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0 

CTi  j 

l  0  / 

The  target  set,  E f,  is  defined  in  terms  of  functional  inequalities, 


(21) 


Ef  '■=  {(X/.C)  e  X  R  :  ef  <  ef(xf,tf)  <  elj] 


(22) 


12  OF  36 


where 


e/(x/>^) 


zf 

Mf 

Vf  sin  7/ 


(23) 


and  and  are  lower  and  upper  bounds  on  the  values  of  the  function  e/,  given  by 


500  ft 

'  500  ft  ' 

d  = 

0.147 

4  = 

0.153 

—  1500-^ 

L  mm  J 

500-^ 

L  mm  J 

(24) 


This  target  set  is  representative  of  desired  conditions  on  approach  to  landing.  In  this  example, 
the  final  angles  of  attack  and  bank  are  not  specified.  Should  the  vehicle  arrive  at  the  target 
set  in  an  unrecoverable  position  (high  angles  of  attack  and  bank  along  with  large  vertical 
speed)  these  variables  as  well  as  their  associated  acceptable  range  of  values  at  the  end  point 
can  easily  be  incorporated  in  the  definition  of  E/. 

Path  constraints  are  derived  primarily  from  the  design  and  mission  considerations  of  the 
vehicle  and  can  be  written  as 

"  —3  g 

— oo 


< 

nz 

< 

6  9 

Jd 

0 

(25) 


where  Jr,  is  defined  in  Eq.  14.  Equation  25  places  a  constraint  on  normal  force  and  addition¬ 
ally  constrains  the  vehicle’s  path  to  reside  in  a  location  in  Mach- a:  space  where  the  vehicle 
can  be  rotationally  trimmed. 

The  state  constraints  are  given  by 
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It  is  now  apparent  that  once  the  problem  has  been  formulated,  it  can  be  posed  as  a 
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constrained,  nonlinear  optimal  control  problem  of  the  following  type: 


Minimize 
Subject  to 


^[x(-)>u(-)^o  ,tf\ 

±(t)  =  f(x(t),u(t),f) 
u (t)  E  U (t,  x(f)) 
x(t)  e  x(t) 

(x(t0),x(tf),t0,tf)  e  E 


where  U  is  written  in  a  form  to  indicate  state-dependent  control  constraints.  Since  solving 
optimal  control  problems,  particularly  state-constrained  problems,  are  widely  considered  to 
be  difficult,26,27  it  is  no  surprise  that  these  problems  belong  to  the  class  of  difficult  problems. 
Nonetheless,  over  the  last  decade  there  has  been  substantial  progress  in  solving  optimal  con¬ 
trol  problems.  Surprisingly,  much  of  this  progress  can  be  attributed  to  a  new  understanding 
of  the  fundamental  issues  of  solving  optimal  control  problems. 27-32  In  order  to  solve  the 
optimal  control  problem,  the  software  package  DIDO  is  utilized.33 

One  last  point  concerning  the  trajectories  computed  in  this  work  is  that  these  trajectories 
are  open-loop.  The  idea  behind  this  work  is  to  be  able  to  incorporate  the  6  DOF  effects 
of  failures  online  to  compute  new  feasible  trajectories.  This  work  provides  a  step  towards 
solving  the  problem  of  including  6  DOF  effects  in  trajectory  replanning  algorithms.  Even 
though  closed- loop  trajectories  computed  online  has  not  been  the  focus  of  this  effort,  it  is 
envisioned  that  the  work  presented  here  could  be  modified  for  closed-loop  guidance. 


Results 

The  following  section  discusses  solutions  to  the  maximum  downrange  problem  under 
nominal  conditions  and  in  the  presence  of  a  control  failure.  As  discussed  in  the  previous 
sections,  we  analyze  a  control  failure  condition  where  the  bodyflaps  become  locked  at  26°, 
which  is  a  case  that  alters  the  control  surface  induced  lift  and  drag  and  also  places  significant 
constraints  on  the  flight  envelope.  The  trajectories  for  both  the  failed  and  nominal  vehicles 
start  at  the  initial  conditions  described  in  Eq.  21  and  satisfy  the  constraints  of  Eq.  25.  The 
initial  conditions  correspond  to  a  point  in  the  flight  envelope  where  the  vehicle  makes  use  of 
its  aerodynamic  control  surfaces  to  control  the  trajectory.  The  vehicle  experiencing  a  control 
surface  failure  must  observe  the  trim  deficiency  path  constraint,  while  the  nominal  vehicle 
is  unconstrained  by  trim  deficiency  because,  as  shown  in  Fig.  2,  the  vehicle  has  adequate 
control  power  to  trim  over  the  entire  flight  envelope. 

The  downrange  profiles  for  both  trajectories  are  shown  in  Fig.  7.  As  expected,  the  down- 
range  for  the  nominal  vehicle  is  much  larger  than  the  downrange  for  the  failed  vehicle.  For 
the  bodyflaps  locked  at  26°,  a  severe  drag  penalty  is  imposed,  hence,  limiting  the  downrange. 
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Figures  8  and  9  display  the  altitude  and  velocity  profiles  for  the  nominal  and  failed  vehicles 
over  the  entire  flight,  while  Figures  10  and  11  display  the  altitude  and  velocity  profiles  for 
only  the  end  of  the  mission.  The  end  mission  requirements  of  500  ft  altitude  and  Mach 
between  0.147  and  0.153  are  achieved.  The  final  velocity  is  about  167.2—  and  the  final 
altitude  is  500  ft.  The  last  endpoint  condition,  from  Eq.  24,  to  be  met  is  the  sink  rate, 
which  is  a  function  of  flight  path  angle.  Figure  12  shows  the  flight  path  angle  profiles  over 
the  entire  simulation  run,  while  Fig.  13  displays  the  flight  path  angles  near  the  end  of  the 
run.  The  flight  path  angles  for  these  two  trajectories  differ  as  a  consequence  of  the  controls 
available  to  the  individual  vehicles.  The  failed  vehicle  exhibits  a  very  steep  descent  but  ends 
the  trajectory  with  a  similar  terminal  maneuver  as  the  nominal  trajectory.  The  similarity 
of  the  terminal  maneuvers  is  a  direct  result  of  the  endpoint  conditions  shared  by  the  two 
trajectories.  Computing  the  sink  rates  yield  values  of  about  —9^  for  the  nominal  vehicle 
and  3-^  for  the  failed  vehicle.  Therefore,  all  desired  endpoint  conditions  have  been  met  and, 
while  the  time  histories  for  all  states  are  not  shown,  the  state  constraints  are  not  violated. 

The  failed  bodyflaps  do  produce  significant  aerodynamic  changes.  Physically,  the  nominal 
trajectory  trades  excess  speed  for  altitude,  as  can  be  seen  in  Figs.  8  and  9.  At  the  top  of 
the  zoom  climb,  the  vehicle  assumes  the  best  lift /drag  configuration.  The  trajectory  for  the 
vehicle  with  the  locked  bodyflaps  is  prevented  from  flying  the  same  angle  of  attack  profile  as 
the  nominal  vehicle  because  of  the  trim  deficiency  constraint.  The  trajectory  does  use  the 
largest  allowable  angle  of  attack  but  this  is  insufficient  to  produce  the  zoom  climb  observed 
in  the  nominal  trajectory.  This  is  more  evident  when  comparing  the  angle  of  attack  profiles, 
as  shown  in  Fig.  14.  The  nominal  constraints  on  angle  of  attack  of  —10°  to  50°  are  clearly 
not  violated  and  significant  differences  between  the  nominal  and  failed  profiles  exist.  Notice 
that  the  angles  of  attack  of  the  failed  vehicle  are  small,  due  to  the  pitch  deficiency  constraint. 
When  the  speed  of  the  failed  vehicle  is  below  about  Mach  =  0.4,  the  trim  deficiency  goes  to 
zero  for  all  angles  of  attack.  Hence,  even  the  failed  vehicle  can  perform  the  flare  maneuver 
shown  in  Fig.  14.  Since  this  example  is  a  maximum  downrange  case,  the  bank  angle  is 
approximately  zero  and  within  the  ±1°  constraint  as  shown  in  Fig.  15. 

A  novel  aspect  of  this  work  is  the  inclusion  of  state-dependent  constraints  into  the  opti¬ 
mization  problem,  in  particular,  the  trim  deficiency  maps.  Recall  that  it  is  desired  that  the 
reshaped  trajectory  he  in  a  region  in  the  Mach- a  space  where  the  vehicle  can  be  rotationally 
trimmed,  if  physically  possible.  The  restricted  region  in  Mach- a-  space  that  results  from 
consideration  of  moment  trim  deficiency  for  the  vehicle  with  locked  flaps,  imposes  severe 
trajectory  performance  penalties  as  shown  in  Figs.  3  and  4.  Figure  16  shows  the  differences 
between  the  two  trajectories  in  Mach-a  space  as  well  as  providing  visual  verification  that 
the  trajectory  for  the  vehicle  experiencing  flap  failure  does  not  violate  the  imposed  trim 
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deficiency  constraint  of  0.003  and  instead  closely  follows  the  contour  in  regions  where  the 
constraint  is  active.  Finally,  we  note  that  there  are  a  number  of  other  issues  in  fault-tolerant 
trajectory  optimization  that  are  beyond  the  scope  of  the  present  paper;  some  of  these  issues 
are  discussed  by  Shaffer.15 


Conclusions 

By  integrating  recent  advances  in  solving  optimal  control  problems  with  inner-loop  con¬ 
trol  allocation,  it  has  been  shown  that  it  is  possible  to  generate  new  feasible  trajectories  for 
an  air  vehicle  that  has  experienced  control  effector  failures.  Methods  were  developed  that  al¬ 
low  the  effects  of  locked  and  floating  aerodynamic  control  surface  failures  to  be  incorporated 
into  3  DOF  models  used  for  trajectory  optimization.  Once  a  failure  has  been  identified,  its 
effects  on  the  ability  of  the  vehicle  to  trim  can  be  determined  and  a  trim  deficiency  map  can 
be  constructed.  For  the  cases  investigated,  the  map  is  a  function  of  Mach  number  and  angle 
of  attack.  Since  Mach  number  is  a  state  variable  (i.e.  scaled  version  of  speed),  incorporating 
trim  deficiency  into  the  optimization  problem  generates  a  state-dependent  control  constraint. 
Although  such  problems  are  notoriously  difficult  to  solve,  results  show  that  pseudospectral 
methods  are  capable  of  solving  them.  When  the  methods  were  applied  to  an  example,  it 
was  shown  that  a  locked  control  effector  can  have  a  significant  impact  upon  the  selection 
of  an  optimal  trajectory  and  that  attempting  to  fly  a  trajectory  designed  for  the  nominal 
vehicle  would  not  be  feasible.  It  was  shown  that  by  recomputing  the  trajectory  using  an 
aerodynamic  database  that  accounted  for  control  surface  failure  effects  and  new  flight  en¬ 
velope  constraints  that  arise  as  a  result  of  the  control  surface  failure,  one  can  successfully 
generate  trajectories  that  can  enable  vehicle  recovery. 
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Figure  1:  X-33. 
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Figure  2:  Pitch  Deficiency  For  Nominal  Vehicle. 
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Figure  3:  Pitch  Deficiency  For  Failed  Body  Flaps  at  26°. 
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Figure  4:  Three-Dimensional  Pitch  Deficiency  For  Failed  Body  Flaps  at  26°. 
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Figure  5:  Nominal  and  Bodyflap  Failure  Lift  Map. 
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Figure  6:  Nominal  and  Bodyflap  Failure  Drag  Map. 
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Figure  7:  Downrange  Profile. 
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Figure  8:  Altitude  Profile. 
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Figure  9:  Velocity  Profile. 
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Figure  11:  Velocity  Profile. 
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Figure  12:  Flight  Path  Angle  Profile. 
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Figure  13:  Flight  Path  Angle  Profile. 
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Figure  15:  Bank  Angle  Profile. 
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Figure  16: 


Trajectory  Trim  Deficiency 


2 


Trim  Deficiency  Map  With  Max  Downrange  Trajectories. 


36  OF  36 


